Metal-insulator transition in the double perovskites 
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We construct an effective Hamiitonian for the motion of electrons among the transition metal ions 
of ordered double perovskites like Sr2FeMoC>6, in which strong intra-atomic Coulomb repulsion U 
is present in only one of the inequivalent transition metal sites. Using a slave-boson formalism, 
we construct a phase diagram which describes a charge transfer transition between insulating and 
metallic behavior as the parameters of the model are changed. The parameters for Sr2FeMo06 are 
estimated from first-principles calculations and a transition to the insulating state with negative 
pressure is obtained. 
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Two mechanisms contribute to the magnetoresistance 
(MR) of Mn Perovskites: one of them is an intrinsic 
mechanism that results from the quenching of spin scat- 
tering of the carriers by localized spins. It dominates the 
high field MR and is most effective at temperatures of 
the order of the Curie temperature T c . The other mech- 
anism which is most effective at low fields, is due to the 
lowering of spin scattering at interphases between two 
regions in the material that have different orientation of 
their magnetization. This mechanism is also temperature 
dependent and is larger at temperatures well bellow T c , 
where the polarization of carriers is large. This explains 
the interest on materials where T c is appreciably larger 
than room temperature. 

The report by Kobayashi et al. |l[ that the double per- 
ovskite S^FeMoOg with T c of about 450 K is half metal- 
lic, with an appreciable low field magnetoresistance has 
renewed interest in these compounds. Furthermore, they 
open new questions concerning the nature of the elec- 
tronic structure. 

In an ionic picture, assuming that all O ions are 2 ~, 
one can imagine the 3d transition metal ions as triply 
or doubly ionized, and Mo to be Mo 5+ (4d 1 ,5" = 1/2) 
or Mo 6 +{4d°,S = 0). In this picture Fe 3+ (3d 5 ,S = 
5/2), would order antiferromagnetically with Mo to pro- 
duce a magnetization saturation of 4fis as observed or 
Fe 2+ (3d 6 ,S = 2), could order ferromagnetically to pro- 
duce the same magnetization. In the extreme situation 
with no carriers on the Mo ion (Fe 2+ and Mo 6+ ), one 
would expect that the system is insulating, because of 
the strong intra-atomic Coulomb repulsion in Fe (U ~ 7 
cV H) compared with the effective Mo-Fe hopping (V ~ 
0.39 eV from our fits explained below or V ~ 0.25 cV ||). 
Actually, there is a certain degree of covalency between O 
and the transition metal ions || . Applying perturbation 
theory to this insulating state, the Mo valence becomes 
vmo = 6 - 24[t Afo _ /(eM - eo)] 2 , where t tm -o is the 
hopping between transition- metal t2 g and O p n orbitals, 
H and ex is the on-site energy for atom X. The mag- 
netic moment of Mo \iMo remains zero in second order in 



tMo-o- This shows that /Umo> rather than vmo can indi- 
cate if the system is insulating or near a metal-insulator 
transition. In fact, the observed Mo-3d chemical shift is 
practically identical to that of M0O3, but vuo < 6 [0. 
For Fe, v Fe - 2 = 4 - fi Fe = 8{t F e-o/ (^Fe - £o)] 2 - 

While Sr2FeMo06 is in fact metallic, neutron diffrac- 
tion experiments obtained magnetic moments consis- 
tent with an insulating state: [iMo = ± 0.1 [Ib and 
fip e = 4±0.1 fiB, indicating that the system is near a 
metal-insulator transition. Our band structure calcula- 
tions give a metallic state with \xmo — 0.2 [1b, and with 
opposite sign as /iFe- The fact that the experimental 
fiMo is lower, points out that the real system is nearer 
to the insulating state than the predictions of these cal- 
culations. Furthermore, substitution of Mo by Re seem 
to increase T c but renders the material insulating, py . 
In addition, substitution of Fe by Co or Mn for example 
makes the compounds antiferromagnetic and insulating 
m. Also, the physical properties (conductivity, magneti- 
zation) are very sensitive to sample preparation. Thus, 
it is of great interest to understand when we can expect 
a metallic or insulating behavior in similar systems. 

In this paper, we intend to address the general question 
of when one should expect insulating or metallic behavior 
in the ordered double perovskites structures. 

In order to investigate the physical properties of these 
materials it is necessary to gain knowledge in their inti- 
mate electronic structure, including the effects of correla- 
tions. We build a tight binding Hamiitonian to describe 
their electronic structure. This Hamiitonian is based on 
the calculated energy bands, as explained later, and is re- 
duced to the minimum set of relevant parameters. Since 
the low-energy properties of the compounds are deter- 
mined by the bands crossing the Fermi energy, the only 
interesting orbitals are the ti g orbitals in the interpen- 
etrating simple cubic lattices of Fe and Mo. Itinerant 
electrons tunnel from the three ti g orbitals of Fe to the 
same ti g orbitals of nearest neighboring Mo, and between 
Mo and Mo with transfer energies V and V'. In the fer- 
romagnetic phase, which will be studied here, only spin 
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down electrons can jump since the spin up orbitals at 
each Fe site are already filled. We consider these orbitals 
to be frozen. The same approach has been followed in 
a recent paper ||. The resulting effective Hamiltonian, 
which contains strong correlations, is treated in a slave- 
boson approximation P,^p[. 

For our band structure calculations we use the Full- 
Potential Linearized Augmented Plane Wave method 



(FP-LAP W) 1 1 1 1 . In brief, this is an implementation of 
density functional theory using the generalized gradient 
approximation (GGA) of the exchange and correlation 
potential. |l3| The Kohn-Sham equations are solved us- 
ing a basis of linearized augmented plane waves. [fL2| Lo- 
cal orbital extensions to the LAPW basis JIJ] are used to 
describe the 3s and 3p orbitals of Fe, the 4s orbitals of Mo 
and Sr. We use a well converged basis set of around 1260 
plane waves and a sampling of the Brillouin zone (BZ) of 
343 points, corresponding to 20 in the irreducible wedge 
(IBZ). We use a muffin-tin radius of 2.2 Bohr for Fe, 2.0 
Bohr for Mo, 2.3 Bohr for Sr, Mo and 1.5 Bohr for O. 

To construct an effective Hamiltonian for the descrip- 
tion of the motion of electrons among transition-metal 
ions, we consider a simple cubic lattice occupied by these 
ions, of lattice parameter a. This lattice consists of two 
interpenetrating f.c.c. sublattices, one occupied by the 
Mo ions, and the other by the Fe ions, in such a way 
that the nearest neighbors (NN) of each Mo ion lie in the 
other sublattice and vice versa. This may be extended to 
other transition metals. In between each two transition- 
metal ions lies an O ion. The O degrees of freedom can be 
eliminated through a low-energy reduction process (see 
for example, Ref. ||[l5| and references therein). If the dif- 
ference between on-site energies tMo — eo and (-Fe — eo, 
is large compared with the absolute value of tuo-o and 
tp e -Oi this can be done by perturbation theory. The xy 
orbitals acquire an effective hopping V to the xy orbitals 
of the four NN in the x — y plane due to a process of 
second order in the hopping between xy and p^ orbitals 
[@|. There are two third-order processes which involve 
0-0 hopping between two p^ orbitals, which lead to a 
next-NN hopping V' between Mo sites. The band cal- 
culations suggest that the equivalent process between Fe 
ions is not important and we neglect it. By symmetry, it 
is clear that the xy orbitals cannot hop to NN or next-NN 
yz or zx orbitals. In addition (excluding direct hopping 
involving d orbitals at distances larger than a/2), the xy 
orbitals at any site i cannot hop out of the plane, because 
its hopping to any O orbital of the sites i± (a/2)z is also 
zero by symmetry. Thus, with a high degree of accuracy, 
electrons occupying xy orbitals move only in the x — y 
plane. Similar considerations extend to the yz and zx 
orbitals. 

This leads to the following effective Hamiltonian for the 
movement of electrons with spin down among d orbitals: 



H = E Fe ^ n ia + E M o ^ n 3 a + U X! 
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Here c ia creates an electron at the orbital a = xy, yz or 
zx of site i with spin down, the sum over i (j) extends 
over the Fe (Mo) sites, rii a — c\ Q Ci a , and S a (j a ) are 
vectors which connect a site with their four NN (next- 
NN) sites lying in the a plane. U is the on-site Coulomb 
repulsion at the Fe sites. The corresponding term at 
the Mo sites is neglected. Eliminating the d — p^ hop- 
ping through a canonical transformation, the values of 
the other parameters are: 

At 2 At 2 

4l Fe _ ^Mo-O 
H>Fe — ^Fe , &Mo — (-Mo — 



tFe-OtMo-O ( 1 
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where to-o is the absolute value of the hopping between 
the p^ orbitals of two nearest-neighbor O atoms. 

In order to treat the Hamiltonian Eq. we use a 
simple extension of the slave-boson theory of Kotliar and 
Ruckenstein ||, to the case in which there are two sites 
per elementary cell and three "colors" per Fe site (instead 
of two representing spin up and down). The Fock space 
at each Fe site i is enlarged introducing boson states rep- 
resented by the creation operators ej (empty), Si a (singly 
occupied at orbital a), d ia p = d^ a (doubly occupied at 
orbitals a and /3 ^ a) and t t (triply occupied). In the 
combined space H reads: 

+ 3*0 

io. ja i a</3 

ia5 Q j a 1a 

+ ^ Ka{s\ a S la + C 4a/3 ( W + A.U - n ia) 

+ E X 'i( e i e i + AU + E 4*Sia + d \ a p d i*p - 1), (3) 
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where z ia represent creation of an electron in the bosonic 
part of the Fock space: 

4 a = (! - e Ui - E " 4yr,<W~ 1/2 

X ( s L e * + E d \ a p S iP + 4 d i-/v) 

x (1 - s\ a s ia - ]T d\ ali d ia0 - tlu)- 1 ' 2 , (4) 
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where 7 and r\ are the two orbitals different from a, anc 
the last two terms of Eq. (||) were introduced to satisfy 
the constraints of vanishing of the corresponding expres- 
sions between brackets, when the energy is minimizec 
with respect to the Lagrange multipliers Xi a and A^ . Th< 
roots in Eq. (Q) were carefully chosen to satisfy three re- 
quirements: i) they are equal to 1 when treated exactlj 
inside this expression, ii) they respect electron-hole sym- 
metry, and iii) when {7 = 0, the saddle point of Eq. (||~ 
reproduces the exact ground state energy. 

In the following we consider the Hamiltonian (|^) ir 
the saddle point approximation, looking at the homoge- 
neous symmetric solution of the minimization problem 
in which all values of the condensed bosons are indepen- 
dent of site and orbital. The problem reduces to a nor 
interacting one with renormalizcd NN hopping Vz, wher* 
z is the saddle point value of z\ a (assumed real). We ne- 
glect the term in V'. Its effect near the metal-insulato] 
transition is small and can be taken into account renor- 
malizing V. Taking V' = allows us to obtain analytica 
expressions for the electronic spectral densities of states 
p a Fe and p a Mo as a function of the saddle point value; 
of the bosons (e, s, d and t) and multiplier Xi a = A: 

PaFe(oj) = {U3- E M o)F(u>), p aM o{u) = (w - E' Fe )F(oj), 



with E't 



Fe 



Ef 



A, F(u) = po(r)/r, 



r = sgn(2co - E Mo - E' Fe )^(uj - E Mo )(uj - E' Fe ), 

and po(u>) is the spectral density for a square lattice witl 
NN hopping Vz. For simplicity, we replace this densitj 
by a constant p = l/(2W) extending from — W to W 
with W — AVz. This allows to perform the integrals 
analytically and the energy per Fe site takes the form: 



E(s,d,t) 



E Mo + E Fe 3[W 2 + (E Mo -E Fe ) 2 /4} 1 / 2 
2 2 

+ [4W 2 /9 + (E Mo - £k) 2 /4] 1/2 

- 6m){E Mo - E Fe - A) + 3U{d 2 + t 2 ), 

where e is eliminated using the constraint 1 = e 2 + t 2 + 
3(s 2 + d 2 ), m = s 2 + 2d 2 + t 2 , and A is eliminated from: 
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Mo 



E> 



AW 
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[W 2 + (E Mo -E' Fe ) 2 /A] 1 / 2 + W \ 
[AW 2 /9 + (E AIo - E' Fe Y/A]^ + 2W/3J 



(7) 



In Fig. 1 we show the resulting values of z after mini- 
mization of E(s, d, t), as a function of U for several values 
of A = Emo — E Fe . For sufficiently large values of these 
quantities, z — 0. This indicates that the carriers be- 
come extremely heavy, the effective band width is zero, 
and the system is insulating. 
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FIG. 1. Reduction factor of the effective hopping as a func- 
tion of U for several values of A = Emo — E Fe . 

For z — > (W — > 0), the expressions (||) and (0) can be 
expanded and the minimization problem can be further 
simplified, allowing a simple analytical description of the 
metal-insulator boundary. Differentiating the resulting 
E(s, d, t) leads to t — and: 
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where y — d/(l — 3m) 1 / 2 . As y varies from zero to 00, 
Eqs. (||) map the metal- insulator boundary. Due to the 
fact that t goes more rapidly to zero than d and m— 1/3, 
the problem takes a similar form as the metal-insulator 
transition in the cuprates |l0| . Taking the limits of large 
or small y in Eqs. (g), asymptotic analytical expressions 
for the boundary can be derived: 
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(9) 



One can see from here that the insulating behavior 
due to Coulomb repulsion at Fe sites is not possible if 
E Mo ~ E Fe < (4^58/3)^ ~ 8.22^, slightly larger than 
the band width of po- This explains the different behav- 
ior of the curves in Fig. 1, depending on the value of 
(E Mo - E Fe )/V. 
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FIG. 2. Phase diagram of the model in the Emo — Efc vs U 
plane, separating the regions of metallic (M), and insulating 
(I) behavior. The dashed lines correspond to the asymptotic 
analytical expressions (^) for U — > and U — > oo. The circles 
correspond to the points listed in Table I. Open circle cor- 
responds to the experimental system at zero pressure. The 
arrow denotes the direction of increasing pressure. 

The resulting metal-insulator phase diagram and the 
asymptotic expressions @ are represented in Fig. 2. The 
realization of such a phase transition depends on the pos- 
sibility of controlling the parameters of the Hamiltonian 
in these materials. As an example, we have studied the 
effects of pressure (positive or negative) on the parame- 
ters for Sr 2 FeMoC>6. The Madelung potential at Fe sites 
should be very different of that at Mo sites due to the 
large charge difference. In consequence, one expects the 
application of pressure or chemical pressure should mod- 
ify the energy difference between Fe and Mo (EMo — Ep e ) 
, on one side and also the hybridization parameter V . 
These parameters were obtained from a fit of our calcu- 
lated t2g energy bands for different values of the lattice 
parameter, to the bands that result from the Hamilto- 
nian Eq. (]l|) with V' = 0, treated in the Hartree-Fock 
approximation. In particular the band that crosses the 
Fermi energy was adjusted at its bottom (point T), and 
Emo ~ Ep e was determined from a fit of the bands at 
the L point. The fit for the experimental lattice param- 
eter is shown in Fig. 3. The bands at higher energies 
are affected by high energy states not included in Eq. 
(||), and a fitting of them is outside our scope. In Table 
I, we show the resulting parameters, and the pressure p 
(obtained from the numerical derivative of the total en- 
ergy in the band structure calculation with respect to 
the volume) for different lattice parameters a. For the 
experimental a = 14.91 A at p = 0, we obtain a small 
but non-zero p due to the errors of the method. The cor- 
responding points in parameter space are represented by 
circles in Fig. 2. One can see that in fact S^FeMoOg 
is not far from a metal-insulator transition and this can 



be achieved applying a negative pressure estimated in -8 
GPa. 
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FIG. 3. Calculated tig energy bands (full line and solid cir- 
cles) and fitting using Eq. (Q) with V' — in the Hartree-Fock 
approximation (dashed lines). The Fermi energy is at 
eV. The wave vectors shown are: W= (w/a, 0, 2n/a), 
L= (ir/a,ir/a,ir/a), T = (0,0,0), X= (0,0,27r/a), and 
K= (37r/2a,37r/2a,0). 



o(A) 


P(GPa) 


V(eV) 


E Mo -E Fe (eV) 


16.51 


-12.47 


0.2417 


2.967 


15.71 


-8.37 


0.3015 


2.886 


15.31 


-4.57 


0.3396 


2.824 


15. If 


-1.69 


0.3618 


2.789 


14.91 


1.72 


0.3850 


2.727 


14.81 


3.52 


0.3964 


2.681 


14.71 


5.36 


0.4096 


2.577 


14.61 


7.17 


0.4203 


2.385 


14.51 


9.07 


0.4494 


2.029 



Table I. Pressure and parameters of Eq. 
ent lattice parameters. 
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In summary, combining band structure calcula- 
tions, with a slave boson technique (equivalent to the 
Gutzwiller approximation) || to treat strong correla- 
tions, we have studied the possibility of a metal-insulator 
transition in Sr 2 FeMoOg.We expect that this approach 
can be used to understand the metallic or insulating be- 
havior in other double perovskites. 
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